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We show that the chameleon scalar field can drive the current phase of cosmic acceleration for 
a large class of scalar potentials that are also consistent with local tests of gravity. These provide 
explicit realizations of a quintessence model where the quintessence scalar field couples directly to 
baryons and dark matter with gravitational strength. We analyze the cosmological evolution of the 
chameleon field and show the existence of an attractor solution with the chameleon following the 
minimum of its efi'ective potential. For a wide range of initial conditions, spanning many orders 
of magnitude in initial chameleon energy density, the attractor is reached before nucleosynthesis. 
Surprisingly, the range of allowed initial conditions leading to a successful cosmology is wider than 
in normal quintessence. We discuss applications to the cyclic model of the universe and show how 
the chameleon mechanism weakens some of the constraints on cyclic potentials. 



I. INTRODUCTION 

A host of observations concord with the existence of a dark energy component with negative pressure, accounting 
for more than two thirds of the current energy budget. The evidence comes, for instance, from measurements of 
the cosmic microwave background temperature anisotropy Q and Type la supernovae 0]. While the data is so far 
consistent with the dark fluid being a cosmological constant, it is nevertheless interesting to consider the possibility 
that near future observations will reveal that w differs from —1. 

Having w ^1 implies that a parameter of the effective Lagrangian, namely the vacuum energy, is time-dependent. 
It follows from general covariance and locality that it must also be a function of space; in other words, the vacuum 
energy is a field, assumed for simplicity to be a fundamental scalar ^. Scalar field models of dark energy generally 
come under the label of quintessence [3|- Of course, this argument assumes that gravity is described by General 
Relativity (GR) for all relevant scales, and it is conceivable that the observed acceleration could result from a break 
down of GR on large scales . However, we focus on the former possibility. 
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Moreover, since w 7^ —1 today, the vacuum energy must have varied significantly over the last Hubble time H^^. 
This in turn requires (f> to have a tiny mass of order Hq ^ 10^^'^ eV. Indeed, if the mass is much smaller than Hq, 
then the field evolution is overdamped and the corresponding w is unmeasurably close to -1; similarly, if the mass is 
much larger than Hq, the field is rolling too rapidly to cause cosmic acceleration. A natural question then arises: if 
such a nearly massless field exists, why have we not detected it in local tests of the Equivalence Principle (EP) and 
fifth force searches 0|? It is well-known that, in effective theories from string theory, such scalars generally couple 
directly to matter with gravitational strength, leading to unacceptably large violations of the EP. 

Recently, two of us (JK and AW) have proposed a novel scenario which offers a natural resolution to this confiict. 
In this work they propose a scalar field which can evolve on a Hubble time today and cause cosmic acceleration, while 
coupling to matter with gravitational strength, in harmony with general expectations from string theory. The basic 
idea is that the scalar field acquires a mass which depends on the local background matter density. On Earth, where 
the density is high, the Compton wavelength of the field is sufficiently short to satisfy all existing tests of gravity; in 
the solar system, where the density is many orders of magnitude smaller, the Compton wavelength is larger than the 
size of the solar system; in the cosmos, where the density is tiny, the field can have a mass of order Hq and cause 
cosmic acceleration. Because its physical properties depend sensitively on the environment, such a scalar field was 
dubbed chameleon. While the idea of a density-dependent mass term is not new [loi [lli Il2| , our work is novel in that 
the scalar field can couple directly to baryons with gravitational strength. 

An important feature of the chameleon scenario is that it makes unambiguous and testable predictions for near- 
future tests of gravity in space. This is timely as three satellite experiments (SEE |^^, STEP Q| and GG Q|) are 
in the proposal stage, while a fourth one (MICROSCOPE Q|) will be launched in 2007. In the solar system, the 
chameleon is essentially a free field and thus mediates a long-range force. This force is very weak for large bodies, such 
as the Sun and the planets, therefore leaving planetary orbits nearly unperturbed. This is because of the thin-shell 
effect, detailed in j^. Intuitively, for sufficiently large objects, only a thin shell just beneath the surface contributes 
to the </)-force on a test mass. This breakdown of the superposition principle is a consequence of the non-linear 
self-interactions of (p. 

Typical test masses in the above satellite experiments, however, do not have a thin shell. Therefore, the extra force 
is comparable to their gravitational interaction. The chameleon model hence predicts that MICROSCOPE, STEP 
and GG could measure violations of the EP stronger than currently allowed by laboratory experiments. Furthermore, 
the SEE project could measure an effective Newton's constant different by order unity from that measured on Earth. 
Such outcomes would constitute strong evidence for the existence of chameleons in our Universe. Moreover, it is 
hoped that the real possibility of such surprising results will strengthen the scientific case for these missions. 

In this paper, we discuss the cosmological history of a universe with a chameleon field. We prove the existence of an 
attractor solution, analogous to the tracker solution in quintessence models, which consists of the chameleon following 
the minimum of its effective potential. For a wide range of initial conditions, spanning many orders of magnitude 
in initial energy density for the scalar field, the solution converges to the attractor. While following the attractor, 
the energy density of the scalar field is always subdominant to the matter and radiation, except for when cosmic 
acceleration is triggered. The onset of the acceleration phase depends on the details of the potential. 

We present a wide class of potentials for which acceleration occurs today; that is, for which the chameleon plays the 
role of quintessence. In doing so, we take advantage of the intriguing fact, showed in _£], that the largest value of M 
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FIG. 1: A typical scalar potential for the chameleon field. 



allowed by existing tests of gravity is 10"'^ eV, which coincides with the energy scale of dark energy. The constraint 
of M < 10~'^ eV was derived for the chameleon in 9] completely independently of any cosmological consideration. It 
is therefore remarkable and unexpected that the energy scale of dark energy emerges from a study of local tests of 
gravity. 

Therefore, a natural class of scalar potentials for our purposes are of the form V{(j>) = M^f{(t)/M), with M « 
lO^'^ eV. That is, potentials that involve a single mass parameter M , which we tune to lO^'^ eV, as required by 
observations. We stress that such a tuning is no better nor worse than the usual tuning of the cosmological constant 
or quintessence models. 

The function / must satisfy only two broad requirements. Following |^, we assume that (z) it is of the runa way 
form, and [ii) it diverges at some finite value of (j), which we take to be = without loss of generality. (See [ig 
for an example of a successful chameleon model where neither of these conditions are satisfied.) Secondly, it must 
be flat and of order unity for today's value of the field, ensuring cosmic acceleration now. A fiducial example is 
V{(t>) = exp(M"/0"), with n some positive constant, which diverges at = and tends to Af^ for (j) ^ M. See 
Fig. n Notice that for pure power law potentials running away to zero at infinity, present-day cosmic acceleration 
with M = eV cannot be obtained. Keep in mind, however, that it is not necessary for V to tend to a constant 
as — !■ cx), as is the case for this particular choice of potential. For instance, the potential could become negative for 
larger field values, as illustrated in Fig. 3 Potentials of this form are of interest because of their direct applicability 
to cyclic models of the universe 

For potentials satisfying the above conditions, we show that the chameleon joins the attractor before the present 
epoch for a wide range of initial conditions. The resulting cosmology is then shown to be consistent with current 
observations. The most stringent constraint on initial conditions comes from the time- variation of coupling constants 
and masses since big bang nucleosynthesis (BBN). Our analysis takes into account 'kicks' due to species becoming 
non-relativistic in the radiation era. We find that the chameleon must be at the minimum by the time of BBN or else 
the electron kick would induce an unacceptably large variation of particle masses. This requires that the chameleon 
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varies by less than 0.1 Mpi over the entire history of the universe. Translated in terms of the initial fractional energy 
density in the chameleon, ri^"*, this requires Q^^^ < 0.1, which is easily satisfied if one assumes equipartition at 
reheating. In particular, we see that this restriction on initial conditions is actually less restrictive than in normal 
quintessence |24j . 

While the cosmological evolution described below somewhat overlaps with previous studies of models of interacting 
d.,- e„e,g, .„d da* UWU (henceforth DE-DM .„de.,, .h„ „e sig„*a.. d«e,e„ces. The 

most important difference is that the chameleon not only couples to dark matter but also to baryons. Thus the 
chameleon model is subject to tight constraints from fifth force and EP experiments. These require that the mass 
scale M in the potential be much smaller than generally considered in DE-DM models. This in turn results in weaker 
constraints from the cosmological evolution of the scalar field. Secondly, the coupling between the chameleon and 
matter fields is exponential (as in as opposed to linear (as in Thus, for small values of the chameleon field, 

this leads to weaker constraints from time- variation of masses and coupling. Thirdly, to our knowledge the discussion 
of the attractor solution and the approach to the attractor in terms of overshoot and undershoot solutions is new and 
parallels the corresponding treatment for quintessence . 

In Sec. ^ we review the main ingredients of chameleon cosmology introduced in P| and give a brief review of how 
existing tests of gravity lead to the constraint M < 10^^ eV. In Sec. lIIII we focus on the cosmological evolution of the 
field for the above class of potentials and show the existence of an attractor solution. How the attractor solution is 
approached for general initial conditions is the subject of Sec. IIVI as in usual quintessence, we find two broad classes 
of solutions, so-called "overshoot" and "undershoot" , corresponding respectively to whether the field begins to the left 
or to the right of the minimum of its effective potential. We derive in Sec. the range of initial conditions allowed by 
cosmological constraints and find that it spans many orders of magnitude in initial energy density of the chameleon. 
It turns out that the main constraint comes from the BBN bound on the time variation of particle masses. The 
behavior of the chameleon during infiation is considered in Sec. I VII We investigate in Sec. I VIII whether the chameleon 



can account for a time- varying fine-structure constant that has been suggested recently [28j . In Sec. IVIIII we apply 
our results to the cyclic model of the universe and show how the chameleon mechanism greatly expands the class of 
potentials suitable for cyclic cosmology. 

II. REVIEW OF THE CHAMELEON MODEL 

The action governing the dynamics of the chameleon field 4> is of the general form 



where Mpi = (SttG)^ ' is the reduced Planck mass, g is the determinant of the metric g^^, TZ is the Ricci scalar and 
'0m'' are various matter fields labeled by i. A key ingredient of the model is the conformal coupling of (j> with matter 
particles. More precisely, the excitations of each matter field ■0m'' follow the geodesies of a metric g^'i which is related 
to the Einstein-frame metric g^i, by the conformal rescaling 

^W^^2ft,/Mp,^^^^ (2) 



where Pi are dimensionless constants |29j . In harmony with general expectations from string theory, we assume that 
the Pi's are of order unity and different for each matter species. Varying the action with respect to yields the 
following Klein- Gordon equation 



Mpi 

where Tj^J = {2/ ^^)6Cn./Sg^^. is the stress-energy tensor for the ith form of matter. 

For relativistic degrees of freedom, it is generally assumed that — 0. This is not quite true however as the trace 
receives two corrections which play an important role in the evolution of the chameleon. First, in a cosmological 
context, each time a particle s pec ies becomes non-relativistic, the trace becomes significantly different from zero for 
about one c-fold of expansion [llj, |30j . A second contribution comes from the trace anomaly HQ Until SecCVI 
we will ignore these two effects and neglect the relativistic fluid contribution to Eq. (j^J. 

For non-relativistic matter with density pi, one can make a perfect fluid approximation to obtain g'^^T^u ~ —Pi- 
Defined in this way, however, pi is not conserved in Einstein frame. Instead, it is more convenient to define a matter 
density pi = p,je^'^''^/^^P' which is independent of and conserved in Einstein frame. We therefore obtain 

^V, + Y: -^P.e^^^"'^^ . (4) 

^ JWpi 

I 

The key realization from Eq. Q is that the dynamics of are not governed solely by V{(f)., but 

rather by an effective potential: 

i 

If V is monotonically decreasing and Pi> (or, equivalently, V{(f)) monotonically increasing and Pi < 0), this effective 
potential has a minimum at 0miii, satisfying 

VA^r.^n) + J2^^p.e^^*-^''^^=0. (6) 
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Meanwhile, the mass of smah fluctuations about (f>aiin is 

= Vf4ct>n.in) = V^04<^mi„) + ^ ^p.e'^'^'"'"/*^- . (7) 

i 

(See [3^ for a stability analysis of the model.) 

The self-interaction potential V{(p) is thought to arise from non-perturbative effects and is assumed to involve a 
single mass scale M: 

V{(P) = mV(0/M) , (8) 

where / is a dimensionless function. As in the original chameleon papers we impose that the potential i) satisfy 
the tracker condition: 

r-^>i, (9) 

and ii) diverge at some finite value of 0, denoted by (/) = </)^. See Fig. In Ref. jol it was believed that these were 
necessary for consistency with current tests of gravity. See |l7| . however, for an example of a successful chameleon 
model where neither i) nor ii) holds. Thus the function / satisfies 

/"/ 

r ' 

/ — !■ oo as a; ^ . (10) 

For most of the paper, we shall assume (p^, = without loss of generality. For applications to cyclic cosmology in 
Sec. IVIlil however, we will need to choose 0* 3> Afpi. 

An essential element of the model is the fact that Vos depends explicitly on the matter density pi, as seen in 
Eq. jsj. In particular, this implies that both (finun and m are also functions of pi. As illustrated in Fig.|31 for a general 
potential satisfying Eqs. the mass is in fact an increasing function of (p: the larger the density, the higher the 

mass. Thus, even though the chameleon mediates a composition-dependent fifth force of gravitational strength, it can 
satisfy the constraints from laboratory tests of the EP and fifth force by acquiring a sufficiently large mass locally. 
This is reviewed below. 



A. Constraint from fifth force and EP experiments 

The tightest constraint on the model comes from searches for a fifth force in the laboratory The potential 
energy U associated with fifth force interactions is generally parameterized by a Yukawa form: 

g-r/A 

U{r) = -aGMiA-h , (11) 

r 

where Mi and M2 are the masses of two test particles, r is their separation, a is the st reng th of the interaction, 
and A is the range. Experiments have found no evidence for a fifth force for A >100 pui |3^, assuming a ~ 0{1). 
Thus, in order for the chameleon model to be consistent with this result, we must impose that the range of the 
chameleon-interaction in the atmosphere, JTiatmi be less than 100 pm. Actually, this ignores the fact that fifth force 
experiments are performed in vacuum, where the density is much less than atmospheric density. Carefully taking into 
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FIG. 3: Chameleon eflective potential for large and small matter density p. This illustrates that, as p decreases, the minimum 
shifts to larger values of (j) and the mass of small fluctuations decreases. The solid curve is the sum of the contribution from 
the actual potential V{(f>) (dashed curve) and the contribution due to the matter density (dotted curve). 

account the presence of the vacuum chamber actually results in a somewhat weaker bound of m'^^^]^ < 1 mm. In 
it was shown that, for the inverse power-law potential V{(j>) — M*^" /(ip — 4>*y^ with n and P of order unity, this 
constraint on matm translates into the following bound on AI: 



M < lO"-" eV . 



(12) 



See the Appendix for details. This condition not only ensures that the chameleon is consistent with fifth force searches, 
but it is also sufficient to satisfy all known local tests of GR, such as EP violation searches, Lunar Laser Ranging 
measurements, tests of post- Newtonian gravity and constraints on the spatial variation of coupling constants 

It is remarkable that the upper bound on AI in Eq. H12|l exactly coincides with the energy scale associated with the 
dark energy today. This comes as a complete surprise since the derivation of Eq. H12|l is based on local tests of gravity 
and is completely independent of cosmic acceleration. We will make use of this intriguing coincidence in turning the 
chameleon into a quintessence field driving the current phase of cosmic acceleration. 



B. Fiducial potential 



Unless stated otherwise, we henceforth take 0* = 0. Since we are interested in applications to quintessence, as our 
fiducial potential we choose 

y(0) =M*exp(ArV<^") (13) 

with M = 10-3 eV. In the limit > M, this reduces to V{(f>) « + M*+"/4>". Since the constant term is only 
relevant dynamically on cosmological scales today, it can be dropped for the analysis of the tests of gravity in the 



8 

laboratory and solar system. Thus Eq. H12() derived in P| for the inverse power-law potential is virtually identical for 
our fiducial potential and hence satisfied for M = 10~'^ eV. 

It is illustrative to apply some of the general expressions above to our fiducial potential. For instance, since 
Vcf, — — n0^"^^M"F in this case, the field value at the minimum (Eq. ©) can be re-written in a form that will be 
useful for later applications 

^min 

where we have assumed a single matter component with density pm and coupling /3, for simplicity. Meanwhile, the 
mass of excitations about the minimum (Eq. lO) is given, once again for a single matter component, by 

where we have used Eq. ©. 

The field value at the minimum today satisfies (fi^^J^ ^ Mpi. This is most easily seen by substituting p^e'^'^'"'"/^^^' ^ 
^('/'min) ~ M'^ in Eq. (|14|) . Assuming (3 and n are of order unity, we find 

/ n,T \ "/("+!) 

Furthermore, since (/),nin is an increasing function of time, the inequality (/jmin ^ Mpi holds for all relevant times, from 
the big bang until today. It follows that the mass in Eq. p5|l can be approximated by 



Mpi < A//pi . (16) 



- K..(^n.n) = ^^^^^^^ |1 + . + n ^- j I . (17) 

Finally, it will prove useful to determine the critical matter density pcrit for which (fimm is of order M. That is, for 
P Pcrit, one has (p^in S> M, and the fiducial potential can be approximated as V{(t)) w M'*-|-M^+"/0". Substituting 
0min ~ -^'^ ^ -Wpi and y(0min) ~ M'^ iu Eq. (|14() . and solving for pcrit, we find 



For (3 and n of order unity and M — 10~^ eV, wc find pcrit ~ 10~®^ ^pp corresponding in cosmological terms to 
a temperature of Tcrit ~ 10 MeV and fl^ w 10^^, where f^m = /Om exp(/3(/)/Mpi)/3-ff^Af|[ is the fractional energy 
density in non-relativistic (matter) degrees of freedom. 



III. COSMOLOGICAL EVOLUTION WITH A CHAMELEON FIELD 



In a homogeneous, isotropic and spatially flat universe, described by the Friedmann-Lemaitre-Robertson- Walker 
metric ds"^ = ~dt^ + a?{t)dx^, Eq. Q reduces to 

ii> + m^^-V^--^ef^'^'-^p^, (19) 
A/pi 

where dots represent derivatives with respect to cosmological time t. For the moment we neglect the effect relativistic 
degrees of freedom have on the trace of the stress tensor. This will be relaxed in Sec. II VI The Hubble parameter H 
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is determined as usual by the Friedmann equation 

3i?2M|i = + y(^) + p^^PHM,. ^ ^ (20) 

where and denote respectively the energy density in matter and radiation. Both and p^ are conserved in 
Einstein frame: 

^ = -AH 

Pm 

^ = ~4H. (21) 

Note that we have simplified the analysis by focusing on a single matter component. It is straightforward to allow for 
different matter species having different (3i, a generalization which docs not substantially alter the results described 
below. 



A. Attractor solution 

We next show the existence of an attractor solution for the chameleon, consisting of the field following the minimum 
of the effective potential, (f> = 0min(^)- Suppose that the field is initially at the minimum. An instant later, due to the 
redshifting of the matter density, the effective minimum will have moved to a slightly larger field value. Clearly, the 
characteristic time scale for this evolution is approximately a Hubble time, H^^. Meanwhile, the response time for 
the field is given by m~^, the period of oscillations about the minimum. From Eq. 117|l . the ratio of these two time 
scales is then 

\i m ^ H , then the response time, m^^, of the field is much shorter than the characteristic time H^^ over which 
the effective potential varies. In this case the chameleon adjusts itself and follows the minimum adiabatically as the 
latter evolves to larger field values. If, however, m <C -ff, the response time is much larger than -ff^^, and the field 
cannot follow the minimum. Instead, the chameleon starts to lag behind the minimum. 

Below we show that, for n and (3 of order unity and M = 10^'^ eV, we have H from the big bang until today. 
Thus, if the field is initially at the minimum, it will follow the minimum as the latter evolves with time. Moreover, 
this solution is stable because if the field is slightly perturbed away from the minimum, it will oscillate and quickly 
settle back to the minimum. In other words, the solution — (j)nun{t) is a dynamical attractor. 

It is useful to compare this with usual quintessence, corresponding to setting (3 = 0. Recall that quintessence also 
has an attractor called the tracker solution Q|, (f> = Attack (0: which is defined by the condition 

That is, this corresponds to the field rolling down the potential at such a rate that its mass is always of order the 
Hubble constant. Our chameleon solution, cj) = 0min(Oi however, satisfies m ^ H. Since m is a monotonically 
decreasing function of ip, it follows that 

</'track(0 > '/'inin(i) (24) 
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at any given time t, i.e., at any fixed value of H. In other words, if (3 were zero (usual quintessence), the field would 
be driven to 0track- For non-zero (3, however, the dynamics of are governed by an effective potential with a minimum 
at 0min < 0track, thus preventing the field from reaching (ptmck- 

It remains to show that m ^ H from the big bang (or inflationary reheating) until the present epoch. While this is 
easily done numerically, it is instructive to provide analytical arguments by studying Eq. H22|l in two different regimes: 
(t'<M and (h > M. 



< M: In this hmit, Eq. if^ imphes 



^>^0n^n^. (25) 



During the radiation and matter-dominated era, is a monotonically-increasing function of time. In particular, 
at the time when the universe has Planckian temperature (the worst case scenario for this argument), one has 
rijii w 10^^^. Thus, for all relevant times, we have fl^ > 10^^*. For M ~ lO^'^ eV, Eq. therefore imphes 



,2 



— > 3/3n • 10^ , (26) 
which is much larger than unity for reasonable values of /3 and n. 
• (j>^ M: In this regime, Eq. H22|l reduces to 

^»3«„ + l,^^n„. (27) 

Moreover, since the potential in Eq. H13|) can be approximated as V{(t)) w Af* in this limit, Eq. H14|l can be 



rewritten as 



Vt^ j^. (28) 



Substituting this into Eq. (|27|l . we obtain 



' \ \ M ) M4 / ' ^^^^ 

where we have neglected a prefactor of order unity. Now SH^Mp^ is an increasing function of redshift and, given 
our choice of M — 10^^ eV, is of order M'^ today. Hence, from the big bang until today. 



Recall from the end of Sec. IIIBI that flm ^ 10^® when 0min ^ M. Since ilm is an increasing function of 
time in the radiation and matter-dominated era, it follows that ^ 10~^ in the regime ^ M. Moreover, 
substituting M = lO^^ gy ^ l0-3"Afpi, we obtain 

|! » 10l2(2n-l)/(n+l) ^ (3I) 

which is greater than unity for n > 0(1). This proves that m > H in the regime (p ^ M as well. 
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Thus we have shown that (f> = (t>min{t) is an attractor from the big bang until today. What will happen to the 
chameleon in the future? Currently, since m 3> -ff , the field is still at the minimum of the effective potential, while 
the universe undergoes cosmic acceleration. Very soon, however, the energy density becomes completely dominated by 
the vacuum energy {V « M^), and the cosmic evolution is driven towards de Sitter. In the process, the matter energy 
density redshifts away at an exponential rate. Eventually, it is so dilute that one reaches a point where to ~ , as 
seen from Eq. H22|l . When this happens, the field can no longer follow the minimum. The dynamics of </> are then 
determined to a good approximation solely by its potential, 1^(0), and the evolution of the chameleon converges to 
that of normal quintessence. 



B. Dynamics of along the attractor 

Having shown that = 0min(i) is a dynamical attractor until the present epoch, we next would like to argue that 
the field is slow-rolling as it follows the attractor and derive an expression for its equation of state. 

It is straightforward to show that 0^/2 <C V as the chameleon follows the minimum. Recall from Eq. © that 0min 
satisfies — V0(0inin) = /3pm exp(/30,„in/A/pi)/Afpi. Taking time derivatives on both sides and using the fact, shown in 
Eq. ((n|), that w V^^ > -V^P/Mpi, we find 

■ _ 3HV^ 

(pmin ~ -r^ 

I/, 00 

It therefore follows that 

9 f VI \ 9 i/2 1 



(32) 



2V{(f>„,in) 2 V^^ \ V,4>^V i 2 to2 r ' 



(33) 



where in the last step we have once again used Eq. (|17|) and substituted F = V4,4,V/V'^ (see Eq. 0). If we recall 
from Sec. lIII Al that > iJ^ and from Eq. ® that T > 1, we obtain 



P ■ 

mill 



« 1 , (34) 



2V{<j)n 

which proves that the field is slow-rolling along the attractor. 

Next we derive an expression for the equation of state of (j) as it follows the minimum. Since (/) is a non-minimally 
coupled scalar, however, this is not given by the usual expression w — (0^ — 2V) / (0^ + 2V) . Instead, we must compute 
w directly from the time evolution of = (jP /2 + V, the energy density in 0. Equation (|34|l implies that ^ V ai 
the minimum, and therefore 

^ « -^0„.i„ = -3ffr-\ (35) 

where in the last step we have substituted Eqs. @ and H32|l . In analogy with the usual energy conservation equation, 
we can define an effective equation of state, Wcs, for as 

^ = -3H{1 + w,s) , (36) 

from which we can read off 

w^g = -l + T-K (37) 
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Equation (|37() is the main result of this Section. It holds for a large class of runaway potentials satisfying Eqs. H10(l . 
Indeed, the only extra assumption in deriving this result is that 0min ^ Mp\ for all relevant times, which implies 

~ V4,4,i(t)min)- 

Applying this result to our fiducial potential, V{(p) — exp(M"/0"), only a few steps of algebra are necessary to 
show that 

l+('l + -V%^V- (38) 



Substituting this in Eq. H37() . we find that Wcs ~ for ^min ^ M, while Wcff ~ — 1 for ^min 3> M. Recall from the 
discussion at the end of Sec. Ill Bi that 0,nin ~ M when the universe has temperature of order T ^ 10 MeV. Thus, for 
this choice of potential, we have shown that the chameleon behaves like dust {wcs ~ 0) for T > 10 MeV, and like 
vacuum energy (wcff ~ —1) for T < 10 MeV. 

Coming back to the case of general V^(0), it is instructive to compare Eq. (|37|l with the usual expression for the 
equation of state, Wusuai = (9^ — 2F)/(0^ + 2V). From Eq. we obtain 

V 1 

Since the second term is down by a factor of /m? <C 1 compared to its counterpart in Eq. H37() . it follows that 
'ii'usuai is much closer to —1 than Weff. But, as argued above, it is really w^s that controls the time-evolution of the 
chameleon's energy density. 

We can also compare Eq. I|37() with the corresponding expression for the equation of state of normal quintessence. 
Along the tracker solution, the equation of state wq of the quintessence field is j3| 

- 2(r - 1) 

" i + 2(r-i) ' (^0) 

where wb is the equation for the background perfect fluid {i.e., « 1/3 during the radiation-dominated era; wb ~ 
during the matter-dominated era). In contrast with wq, the chameleon's effective equation of state is independent of 
wb ■ This is because the location of the minimum in chameleon cosmology is nearly independent of the energy density 
in relativistic degrees of freedom, since these have a nearly traceless stress-energy tensor and therefore contribute 
negligibly to the right-hand side of Eq. |j2Jl. Consequently, the evolution of the chameleon along the attractor is 
insensitive to whether or not the matter density dominates over radiation. 



IV. APPROACHING THE ATTRACTOR 



In this Section we describe how the attractor solution is approached for general initial conditions. For simplicity we 
assume the kinetic energy of the field is initially zero, although it is trivial to generalize the discussion to the case of 
non-zero kinetic energy. Moreover, we focus on the case where the energy density in (j) is initially less than the energy 
density in matter and radiation, as expected from equipartition at reheating. More pragmatically, however, we will 
see in Sec. |V] that this is required in order for the cosmology to be consistent with the measured abundance of light 
elements. 

So let us start the field from some arbitrary value 4>i at time ti. The effective potential at that time displays a 
minimum at 0min(ij): with ^min determined by Eq. ij^Jl. If the field is released at the minimum, (pi — (/'min(ii)' then 
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the field just keeps following the minimum until today, as argued in Sec. IIII K\ Similarly, if the field starts very near 
'^min(^i), it osciUates and quickly settles to the minimum. 

More generally we wish to consider the limiting cases where <^ (/'min(^j) and cjji 3> 0min(ii)- it will soon 
become clear, and to make a parallel with the analogous discussion in usual quintessence, the corresponding solutions 
will be referred to as overshoot and undershoot, respectively. 



A. Undershooting 

In this case 3> (/'min(^i), the source term can be neglected, and the equation for (j) reduces to 

?i + 3i/0«-^T;, (41) 

where we have reintroduced the trace of the stress tensor. If the trace were negligibly small for relativistic degrees 
of freedom, then one would have ~ —Pm, where pm is the non-relativistic matter density. But then, during the 
radiation-dominated era, the latter would be utterly negligible compared to the friction term, '6H(f). Thus the field 
would be overdamped and would remain essentially frozen at its initial value 

Fortunately, however, the trace is not always small for a realistic relativistic fluid. As the universe expands and 
cools, each massive particles species successively becomes non-relativistic whenever m ~ T. When this happens, the 
trace becomes non-zero for about one e-fold of expansion, thus driving the field a bit closer to the minimum jsol . To 
see this explicity, we follow 0,0] and note that each component, labeled by i, of the relativistic plasma contributes 



where g,{T) - Ebosons 5^''°"(r./r)4 + (7/8) Efcrmions 5f ™'°"(7^./r)4 is the usual expression for the effective number 
of relativistic degrees of freedom, while gi and Tj are the number of degrees of freedom and temperature of the i^^ 
species, respectively. The function r is defined by 



r{x) ^ f du^^^^ , (43) 

where the ± sign is for fermions and bosons, respectively. It is negligibly small both for a; ^ 1 and a; 3> 1, but 
becomes of order unity when x ^ C'(l)- 

Integrating Eq. H41|) numerically, one finds that the total displacement in (p due to the source term in Eq. (|42|) is 
approximately given by 

(A</,),«-/?^^|^/^|Mp,, (44) 

where the upper and lower numerical coefficients are for fermions and bosons respectively, and gi,{mi) is the effective 
number of relativistic degrees of freedom when T = rrii. This can also be derived analytically by approximating t(x) 



as a (5-function source in Eq. 1)42(1 . as shown in the Appendix fSec. IXIB)l . Thus, each species, whenever it becomes 
non-relativistic, effectively gives a "kick" to the chameleon, driving the latter closer to the minimum of the effective 
potential. 
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Particle 


9i/g* 


t 


12/106.75 


Z 


3/95.25 




6/92.25 


b 


12/86.25 


T 


4/75.75 


C 


12/72.25 


TT 


3/17.25 


A* 


4/14.25 


e 


4/10.75 



TABLE I: List of particle species, in order of decreasing mass threshold, that contribute to kicking the chameleon towards the 
minimum. 

We can calculate the maximal displacement, (A0)tot, by summing over all relevant massive particle species. In 
doing so, we only include particles that become non-relativistic at a much lower temperature scale than that of the 
phase transition through which they acquired a mass. See Table^for a list of the relevant particles and corresponding 
(7i's. We will see below that big bang nucleosynthesis (BBN) constrains the chameleon to be at the minimum by 
the onset of nucleosynthesis. This will be the case if {(f)i — (j)^^-^^^) < |(A(/))tot|, where (t>m^^^ location of the 

minimum at BBN. Hence, in doing the sum, we must neglect the electron contribution since its mass threshold is 
too close to the temperature during BBN (~ 1 MeV) . With this proviso, we can substitute the results of Table ^ in 
Eq. (O and obtain 

(A0)tot ~ -/3Mpi . (45) 

Thus, for f3 ^ 0{1), the cumulative effect of the kicks is to push the field a distance of order AIp\ towards the 
minimum. 

The results of numerical integration are shown in Fig. ^for (j)i — 2 Afpi, where we compare the solution with and 
without the kicks described above. Because of numerical limitations, we were restricted to choosing (pi > Afpi only. 
For the solution without kicks, the chameleon remains frozen at its initial value throughout the radiation-dominated 
epoch. For the solution including kicks, however, it is pushed to smaller field values as various particle species 
successively become non-relativistic. It is seen that the total displacement from z = 10^° to z = 10^ (redshift of BBN) 
is indeed of order ^ Mpi, consistent with Eq. 145|l . 

As we will see later, the particular solution shown in Fig. 0] is ruled out by BBN constraints since the field is 
not at the minimum by the onset of BBN. This is because numerical limitations forced us to choose > Mpi, a 
range of initial values larger than the total displacement given in Eq. H45|l . Thus, in this case, the field begins so far 
away from the minimum that the sequence of kicks is insufficient to bring it in the vicinity of the minimum by the 
onset of BBN. Even though this particular solution is ruled out, a few instructive comments can be made. After the 
electron has become non-relativistic {z « 10^), the chameleon receives no further "kicks" and thus is essentially stuck 
at some particular value. It remains frozen there until matter-radiation equality (z « 10^). At that time, one has 
iJ^Mpj ^ piji « r^, and the driving term in Eq. Ij41(l is then comparable in magnitude with the friction term. Thus 
the field begins to roll towards the minimum, undergoes large anharmonic oscillations, and eventually converges to 
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the minimum. These oscillations will be described in more detail in Sec. IIV CI 

Had the numerics allowed us to probe the regime (j)i < Mpi, the solution would initially look similar to that of 
Fig. 21 with each kick bringing the field closer to the minimum. Eventually, however, one of the kicks pushes the 
field sufficiently close to the minimum. The chameleon then starts oscillating and quickly settles to the minimum. 
Subsequent kicks generate oscillations about the minimum which are rapidly damped by the expansion of the universe. 
In particular, for the case of the electron contribution which kicks in during BBN, these oscillations are sufficiently 
small in amplitude to obey BBN constraints on time-variation of particle masses. This is discussed in more detail in 
the Appendix. 

To summarize, if > /3Mpi, the chameleon is initially so far from the minimum of the effective potential that the 
sequence of kicks is insufficient to push it in the vicinity of the minimum. In particular, the field is still away from 
the minimum by the onset of BBN, and, as we will see in Sec. this solution is ruled out by BBN constraints on 
time-variation of particle masses. If 4>i < /3Mpi, then the kicks eventually push the chameleon near the minimum. At 
that point, the field oscillates and quickly settles down to the minimum. In particular, it is at the minimum by the 
onset of BBN, and this solution is consistent with BBN constraints, as explained in Sec. 

We conclude by noting that another potentially important contribution to T^^ during the radiation-dominated era 
is from the trace anomaly. It was shown in [31| (see also |33|) that the effective equation of state for a plasma of an 
SU(iVc) with coupling g and Nj flavors is given by 



6tt^ y^TT J 2+ lNcNf/{N^ - 1 

For QCD with Nc = 3 and Nj — 6, the above is of order 10^^ for energies above 100 GeV {i.e., in the perturbative 
regime of QCD). This is too small to yield a significant displacement for the chameleon. However, larger gauge groups 
and Nf could result in an important driving term near the unification scale thereby improving on Eq. H45|l . For 
the purpose of this paper, however, we take a conservative approach and restrict ourselves to Standard Model degrees 
of freedom. 



B. Overshooting 

Since (f>i <^ if'min{ti), in this case it is the source term that is negligible in Eq. 1)19(1 . Thus the equation for can 
be approximated by 

+ 3il<^ « -V^ , (47) 

which describes the evolution of a minimally-coupled scalar field. Once again, the dynamics are governed by a friction 
force, 3H(j), and a driving force, —V^,. However, the driving term initially dominates over the friction term in this 
case. To see this, note that (f>i <C 0min(ii) implies V^00(0i) S> mf, where is the mass of small fluctuations about 
4>min{ti). Moreover, we have argued in Sec. IIII A1 that ^ for all relevant times, therefore V^^^{(l)i) ^ i/f, 
where Hi is the initial value of the Hubble parameter. In other words, the field is undcrdamped, and its evolution is 
essentially that of a free field. 

Thus very quickly its energy becomes kinetic-dominated, (jP' 3> V . Since (f>i is much smaller than ^min, the field 
then rolls past the minimum and thus overshoots. It keeps on rolling until its kinetic energy has sufficiently redshifted 
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Solution including kicks 
Solution neglecting kicks 
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le+16 



le+20 
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FIG. 4: Evolution of the chameleon field as a function of 1 + z. In this example the potential is V{(f>) = /(f>'^ , P ~ ^ and 
M = 10 MeV, the latter being larger than the required M — 10~^ eV due to numerical limitations. The initial conditions for 
the field are =2 Mpi and <f>i = at z = 10^" . The dashed line presents the case in which the contribution of particles which 
become non-relativistic is neglected. The solid line includes the contribution of the particles in TableQto T^. Since (j>i ^ /3Afpi 
in this case, the field does not reach the minimum by the onset of BBN {z ~ 10^). As explained in the text, this solution is 
therefore ruled out by BBN constraints on time- variation of particle masses (see Sec.lVl. 

so that the Hubble damping term becomes important. Then the field essentially comes to a halt at some value 
'/'stop > 0miii, which we can estimate as follows. Since the energy in (/) is kinetic-dominated as it rolls, we have 
(j) ^ a^^. Moreover, for the range of initial conditions relevant to this discussion, we may assume the universe is 
radiation-dominated, and therefore a ~ t^^^. Combining these two facts, a little algebra shows that 



where ft^ and are the initial fractional energy density and initial value of the scale factor, respectively. With 
initial condition (j) = (pi sX a = at, this integrates to 




(48) 




(49) 
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which gives (f) — > (/>stop in the Hmit ai, where 



This is precisely what one finds for the overshoot solution in usual quintessence {e.g., see Eq. (11) of which is not 
surprising since the non-minimal coupling of the chameleon is irrelevant in the above derivation. Once the chameleon 
reaches 0stopj the solution is then exactly as in the undershoot case above, with <f>i replaced by 0stop- 

C. Converging to the minimum 

When the field reaches the vicinity of the minimum, it begins oscillating and eventually converges to the minimum. 
In this Section, we study these oscillations and the rate of convergence. To proceed analytically, we assume that the 
linear approximation is valid so the oscillations are harmonic. Strictly speaking, this of course only holds when the 
field is sufficiently close to the minimum. Nevertheless, we will find numerically that the range of validity is actually 
much wider and applies even when the oscillations are apparently anharmonic. Moreover, we neglect the kinetic 
energy from the last "kick" that drove the field in the vicinity of the minimum. 

In the linear approximation, the effective potential is given by 

V,s~\'m'{t){^{t)-^^,^f , (51) 

where we have dropped the vacuum energy term, which is completely irrelevant for the present discussion. Meanwhile 
the total energy density in (j) is well approximated by 

P4, « \{^ - 0min)' + ]^m'{t){m ~ 0„,in)2 . (52) 

Note that both 0niin a-nd m are functions of time, a fact that will be important in the discussion below. Since the 
motion is assumed harmonic, averaging over a few oscillations gives 

P4 « im^^t) < (0(t) - 0,„i„)2 > . (53) 

Since m ^ H, the problem is analogous to a pendulum which is slowly lengthened [s^l . It is well known that in 
this adiabatic approximation the total oscillator number, 

m(t) ' 

is conserved. Substituting Eq. H53|l . it follows that 



AA^^, (54) 



m{t) < m) - <j)^^Y >~ . (55) 

See jSlijl and the Appendix for alternative derivations. 

If m were constant, Eq. H55I) would imply that the energy density in the oscillations redshifts like dust, which is not 
surprising. In the chameleon model, however, this quantity is time-dependent. Recall from Eq. I|17|l that m w ^Jj^^j 
and thus 

m 2 Vm 2 Vi. 
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where in the last step we have substituted Eq. (|32(l . For the inverse power-law potential, V ^ (j) this reduces to 

m 3Hn + 2 



from which we conclude that 



Substituting this in Eq. H55|l . we obtain 



m 2 n + I 



m ^ . (58) 



(57) 



< 



^/-.ninWf >'/'-«"'"/'^"+'^ (59) 



It is worth mentioning that the time-dependence of < (0 — 0min)^ > is independent of M. The parameter M only 
enters implicitly in that 0miii at any fixed time depends on the choice of M. Thus, for fixed (pi, the difference 
< {(j) — (pmin)'^ > is initially larger for smaller M, and therefore it takes longer for to converge to ^min- 

The above results have been checked numerically, as shown in Fig. |31 The dotted line has the slope predicted by 
Eq. ifS^ . We see that the latter is a good fit throughout the oscillatory regime, even when the oscillations appear to 
be anharmonic. 

Finally, we note in passing that Eqs. lf55|) and ifS^ together imply 

P0 ~ a-3(3n+4)/2(„+l) (go) 

For any n > 0, we see that the energy density in the oscillations redshifts faster than radiation. Therefore, the 
chameleon does not suffer from the old moduli problem. 



V. CONSTRAINTS ON INITIAL CONDITIONS 

Next we use the results of the previous section to determine what subset of initial conditions gives a cosmology 
consistent with current observations. As derived in Sec. IIVI fsee Eq. H45ll l. if the chameleon starts at, or overshoots 
to > (3Mp\, then the field will not be at the minimum of the effective potential by the beginning of BBN. We now 
argue that this case is ruled out by BBN constraints on time- variation of particle masses. 

Due to the conformal coupling in Eq. (O, a constant mass scale m*^*) in the matter-frame is related to a (/)-dependent 
mass scale to(0) in Einstein-frame by the rescaling m(0) = e''''*/*^P'm*-*^. Thus variations in lead to variations in 
the various masses: 

Am 



m 



Nucleosynthesis constrains the variation in m from the time of nucleosynthesis until today to be less than 10% or so. 
Since the value of 4> at the minimum today is much less than Afpi, this constrains (/)bbn, the value assumed by (p at 
BBN, to satisfy 

0BBN <0. 1/3-1 Afpi. (62) 



Evidently, if the field starts at, or overshoots to, ip > /SAfpi, then it will not yet be at the minimum by the onset of 
BBN, assuming f3 ^ 0{1). In this case, (/>bbn violates the above bound. 
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FIG. 5: The root mean square quantity <(</> — (^min)^ >i/2 plotted as a function of redshift for the undershoot solution with 
4>i = 10^^ i\'fpi and neglecting kicks. The potential is as in the previous two figures. The dotted hue is the slope predicted by 
Eq. lEU. 



Thus, for the undershoot and overshoot cases, we have the respective constraints (j^i < /3A/pi and 0stop ^ f3Mp\. 
Using Eq. (|5U|) . the latter translates in a bound on the initial fractional energy density in the chameleon: 



(63) 



where we have neglected and assumed f3 ^ C(l). This is the main result of this Section. Remarkably, this is a 
much weaker bound than in normal quintessence, where fi^'' is required by BBN to be less than 10~^. One would 
have expected the opposite since, after all, the chameleon couples directly to matter whereas normal quintessence 
does not. But, as shown above, it is precisely this direct coupling that pushes the field more effectively towards the 
attractor solution, through a sequence of kicks, thereby resulting in a weaker constraint. The above upper bound on 
ri^'' allows for a range of initial scalar field energy density spanning many orders or magnitude and consistent with 
equipartition of energy at reheating. 

Equation H()3|) ensures that the field is at the minimum from the onset of BBN until the present epoch. In the rest 
of the Section, we show that it is also a sufficient condition to satisfy all current cosmological constraints and thus to 
obtain a successful cosmology. 



A. Cosmic microwave background anisotropy 



A potentially important effect of the chameleon on the CMB is the modification of the distance to the last scattering 
surface. Indeed, Eq. (|61|l implies in particular that the electron mass varies by 
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where (j)^^}^ and (t>[^^^^^ are the field values at the minimum today and at recombination, respectively. Such a variation 
in the electron mass modifies the binding energy of hydrogen and thus changes the redshift of recombination, Zioc, by 
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The WMAP experiment constrains z^cc to within 0.1% or so. From Eqs. (|14|l and we find that 
and 0min/Mpi (A//Afpi)"/("+^\ which trivially satisfies the WMAP bound. Thus the distance to last scattering is 
virtually identical to that in normal quintessence. 

The other effects of the chameleon on the CMB are the backreaction of its energy density and its influence on the 
growth of perturbations. The former is trivially satisfied. Indeed, for temperatures less than 10 MeV (z < lO^*'), the 
energy density in </> behaves effectively like a (small) cosmological constant, as shown in Section III B. We now turn 
our attention to the linear evolution of density perturbations. The full non-linear analysis requires numerical work 
which lies beyond the scope of the present paper. 

B. Density perturbations and large scale structure 



We study perturbations in the synchronous gauge, where the perturbed line element has the form 



(66) 



Throughout this Section, we use conformal time r rather than cosmic time. The perturbation equations have been 
discussed in many papers (see, e.g., ^^nd references therein), so we simply state them here for pressureless dust 
and the chameleon field. We follow the literature and write the equations in Einstein frame. 

The evolution equations for the dark matter density contrast, Sc — SipmS^'^^'^''^^) / PmS^'^^'^''^^ , and the divergence of 
the velocity field of the dark matter fluid, Oc — V-v, are respectively given by 
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(67) 



Mpi' 



(68) 



where primes denote derivatives with respect to t, and h = 5^^hij is the trace of hij. For the latter, Einstein's 
equations give 



h" + aHh' = 



Meanwhile, the perturbed Klein-Gordon equation for the chameleon reads 

(<50)" + 2ai/(d»' + + aV^^) 5^ + = 



/3 2 
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(69) 



(70) 



We can safely ignore terms proportional to (j)' since the field is slow-rolling along the attractor, as shown earlier. Then, 
taking the time-derivative of Eq. (|67(l and using Eqs. (|68|l . H69|l and l|70|l . we obtain the following equation for &c'- 



5'' + aH5' =-a^H^ 



2/32 



1 + aW^^/k^ 



Sc, 



(71) 
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where the quantity in brackets can be interpreted as an effective Newton's constant. In particular, the term propor- 
tional to y 00 results from the chameleon-mediated force 371 , which is negligible if the physical length scale of the 
perturbation is much larger than the range of the 0-mediated force, that is, if a/k ^ • In. this case the left hand 

side of Eq. (|71|l is well approximated by 3a^i7^(5c/2 and the dark matter fluctuations grow as in GR. 

If the field is at the minimum, then m? « 3> i?^, as shown earlier. For instance, our fiducial potential with 
n > 1 gives m/H > 10^" at recombination. Thus the chameleon length scale is much less than the Hubble horizon at 
the time of recombination, which therefore implies that the CMB spectrum is unaffected. This has to be contrasted 
with normal quintessence models, for which the mass of the field is on the order of the expansion rate along the 
tracker solution. In this case, fluctuations on scales of the size of the horizon cannot be neglected and leave distinctive 
features in the anisotropy spectrum of the CMB 0| . 

Writing V{4>) ~ M'^ + M^^"' /(/)" , the length scale Acham = ^^^^^^ below which perturbations feel a different Newton's 
constant is given by 

l+n/2 / (h\ ^+"/^ 

10"^ (tt) cm, (72) 



where in the last step we have assumed n ~ 0{\) and substituted M = 10"'^ eV. Density perturbations below that 
scale will initially grow faster because the effective Newton's constant is larger. From Eq. (|71|) the modified growth 
rate is (5 oc r^, where 



^ = -^±y^ +6(1 + 2/32). (73) 

For P ~ 1, for example, one obtains x = 3.8 for the growing mode, compared with x = 2 as predicted by GR. Such 
perturbations thus grow faster and enter the non-linear regime earlier. From Eq. H16|l with M = 10~^ eV and n > 1, 
we find that Acham at the present epoch satisfies 

Acham < 100 pC . (74) 

While this is a rather small length scale, it would be interesting to investigate if this could change substantially the 
details of galaxy formation, and if early star formation (and therefore early reionization) is a natural consequence 
of chameleon cosmology. Another avenue worthy of investigation is whether the chameleon favors the formation of 
supermassive black holes at the center of galaxies. 



VI. THE CHAMELEON DURING INFLATION 



Let us now study the behavior of the chameleon during a period of inflation in the early universe. During inflation, 
the effective potential for the chameleon field is 

Fcff (<^) « exp(M"/^") + pvace^'^*/*^- , (75) 

where the last term arises from the coupling of the chameleon to the inflaton vacuum energy, pvac- The factor of 4 in 
the exponential follows from the fact that 1 — 3w = 4 in this case. Following the analysis leading to Eq. H18|l . wc find 
that this potential has a minimum at some (padn <C M. Since pvac is essentially constant during inflation, so is 0min, 
and the chameleon is stabilized. 
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The mass of the chameleon at the minimum is readily obtained from Eq. (|22(l by letting /? ^ 4/3 and ^vac ~ 1: 

« Uf3n^ ( ^] » . (76) 

M V / 

As the mass is constant during inflation, the chameleon oscillates around the minimum with an average amplitude 
given by 

<{(b~ ^,rnnf >~ a-^ . (77) 

As expected, the chameleon field behaves like a dust component during inflation. Due to the exponential growth of 
the scale factor, this implies that (/) very quickly settles at the minimum of the potential. 

At the end of inflation when the universe reheats, the inflaton decays not only into radiation and matter, but also 
into coherent kinetic energy of the chameleon. However, the production of chameleon quanta is generally suppressed 
because of their large mass. Since the universe becomes radiation-dominated at the end of inflation, the minimum 
of the effective potential suddenly moves to a much larger value. The chameleon is therefore released from a field 
value much smaller than (/)min- Hence, we conclude that the overshoot solution, discussed in Section [iV Bl is the more 
natural outcome of an inflationary phase. 



VII. CHAMELEON AND TIME- VARYING a 



Recent analysis of absorption spectra of quasars have led some to claim that the fine-structure constant oem might 
have evolved by approximately one part in 10^ over the redshift range 0.2 < z < 3.7 j^l- If this turns out to be true, 
then general covariance would imply that ckem can vary both in space and in time, that is, it must be a function of 
a field. A space-time varying fine-structure constant was first posited by Jordan 0| , Teller , Stanyukovich j4l| , 
while the implications for EP violations were studied by Dicke 42] and more recently by Bekenstein |43| . The recent 
evidence from quasar spectra has rejuvenated this idea, resulting in a flurry of papers |4J,|4a|. 

Here we show that the chameleon cannot lead to a time-variation of oem consistent with the recent observational 
claims. Since we have so far assumed that the chameleon couples conformally to matter fields, as seen from Eq. 
and since the Maxwell action is conformally- invariant, at tree-level the chameleon does not lead to a time- varying 
OEM- However, one can easily consider generalizations to the matter action in Eq. such as 

2^*/*^-An(Vi^\5i.^), (78) 

where /3 is a constant. In particular, the Maxwell lagrangian in this case reads g^^e'^'^/^^^'i^^i/F^", corresponding to 

OEM - e-^^/^^- . (79) 

Then, variations in induce variations in qem of magnitude 

AttEM 



/S^. (80) 

ttEM JWPI 

Since the cosmological value of the chameleon, (/)(z), is a decreasing function of density and, therefore, a decreasing 
function of redshift z, we can approximate Eq. H80() by 

a(0) / A/f \ 
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where we used Eq. (|16|l in the last step. For n > 1, the resulting Aq;em/«em is many orders of magnitude too small 
to account for the time- variation advocated by |28i |. 

Taking into account local variations in the chameleon does not help. Indeed, recall that for sufficiently dense 
objects, the local value of the chameleon depends mostly on the local density and is insensitive to the background 
cosmic density. The absorption clouds used in the analysis of |3| can be classified in three populations: weak 
Mgll systems, Lyman limit systems and damped Lyman- alpha systems, with estimated density of 10~^^, 10~^^ and 
10~24 _ io~^^ g/cm^, respectively. The key point is that all of these are comparable with the local galactic density 
of lO^^"* g/cm"^, and thus the chameleon value in these systems should be nearly the same as locally. 



VIII. APPLICATION TO THE CYCLIC UNIVERSE MODEL 



To illustrate the usefulness of the chameleon mechanism in designing cosmological scenarios, in this Section we apply 
our results to the cyclic model of the universe [isi [igi l2o[ l2l[ bslbsf . The cyclic scenario proposes that time did not 
begin at the big bang, as assumed implicitly in standard inflationary cosmology, but rather extends infinitely far in the 
past as well as in the future. Thus our current epoch of expansion is only one out of an infinite number of cycles. The 
cyclic model addresses the homogeneity, flatness and monopole problems of the big bang model, and generates a nearly 
scale-invariant spectrum of density fluctuations, without invoking a period of high-energy accelerated expansion. As 
such, it constitutes the most serious candidate for a viable alternative to the inflationary paradigm. 

Except for very near the big crunch/bang transition, the cyclic model is well-described by a four-dimensional 
effective action of the form given in Eq. (Q: 

|d4^V^|M7^-i(9^)2_y(^)|_|d4a;£^(^«,g«). (82) 

In particular, as in the simplest inflationary models, the scenario consists of a scalar fleld rolling down its potential 
V{(j)). An important distinction with inflationary theory, however, is that here (f> has a higher-dimensional interpreta- 
tion of measuring the distance d between two end-of-the- world branes (more precisely, orbifold fixed planes). In other 
words, is a radion. In this higher-dimensional picture, the big crunch/bang brid ge b etween each cycle corresponds to 
the cataclysmic collision of the branes. The relation is d = L ln[coth(0/Mpi\/6)] [islll^ . and thus the brane collision 
(d — !■ 0) corresponds to the limit 4> — > oo. 
At tree-level, the metrics g^jll are given by Eq. |2Il with couplings 

A = ^ , (83) 

which corresponds to the Kaluza-Klein limit in the higher-dimensional picture. We henceforth assume that this regime 
holds for the relevant range of (p. The potential V{4>) is thought to arise from non-perturbative effects, such as virtual 
exchange of branes in the higher-dimensional theory. A typical cyclic potential, sketched in Fig. 12 is given by 

V{(t)) = A/''e^'^"/('^-'^*'" {l - e-=('/'-'/'cr==3)/Mp, I . ^ (84) 

where we have reintroduced 0^,, the value of (j) for which the potential diverges. The field value t/fcross is where the 
potential vanishes. The positive constant c must satisfy c > 0(10) in order for the spectrum of density perturbations 
to be nearly scale invariant 21]. The function J- {(I)) accounts for the fact that the non-perturbative effects must turn 
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off as (/) — > cx), since the string coupling goes to zero in this hmit. In order for the Kaluza-Klein limit assumed in 
Eq. 183|l to apply for all (f> > (pi,, we impose that 

0* » Mpi ■ (85) 
Currently, the field lies at cf)'-^^ — 0* ^ Af , indicated by a dot in the Figure. Thus, V{(j)o) ~ M^, and this potential 



energy drives the observed acceleration of the universe (Region a)). This phase of cosmic acceleration lasts sufficiently 
long to empty out our observable universe, thereby making it highly homogeneous, isotropic, spatially-fiat and nearly 
vacuous. After a while, the field begins to roll down the potential and reaches Region b). Since V < Q m this region, 
cosmic expansion eventually comes to halt, and the universe enters a phase of contraction. It is during this phase 
that a nearly scale- invariant spectrum of density perturbations is generated from quantum fiuctuations in cj) |l9l ^ | . 
When the field reaches Region c), the function J-{(l)) in Eq. H84I) becomes important and causes V to go to zero 4G|. 
This ensures that the energy density of the universe is dominated by the kinetic energy in (f) as it zooms toward +00, 
as required by the prescription of |20j. The — > cx) limit corresponds to the big crunch/bang transition, at which 
point the universe reheats and becomes filled with thermal matter and radiation. This marks the beginning of the 
hot big bang phase. Meanwhile, 4> bounces back, rushes through Region b) and eventually comes to a stop in Region 
a). It remains essentially frozen there, until the universe is sufficiently cold to allow the vacuum energy in cj) to drive 
cosmic acceleration. The cycle then repeats itself. 

It is crucial that does not result in EP violations stronger than allowed by experiments. One way to ensure this, 
as proposed by Steinhardt and Turok ^isl, is if the couplings f3i in Eq. 1)83(1 are functions of and satisfy 

A;(0W)<lO-\ (86) 

for today's value of the field, (j)^^'' . This possibility is certainly allowed by some brane- world models, such as the 
Randall-Sundrum scenario 47J. In general, however, one expects that Pi will be of order unity and different for 
different matter species. 

Here we argue that no such extra condition is necessary since (j) is in fact a chameleon. Indeed, the reasoning that 
lead to Eq. 1)86(1 neglected the effect of the background matter density. In Region a) , the potential is approximately 
given by 

V{<j)) w M'*e*^"/('^-'^*'" , (87) 

which is of the same form as our fiducial potential (see Eq. I(13|l '). Moreover, since the various /3i's are non-zero and 
positive, the dynamics of are governed by the effective potential in Eq. © with Pi — 1/V6. 

The above analysis therefore greatly expands the range of models and brane-world set-ups suitable for cyclic 
cosmology. The only requirements are that Pi > and that the potential be of the general form 

where / satisfies Eqs. ((10(1 . In particular, the chameleon mechanism disposes of condition 1(86(1 . 



IX. CONCLUSIONS 



We have explored the complete cosmological evolution of the chameleon field for general models where the effective 
potential displays a minimum, as in the original scenario. Our analysis shows that the chameleon can act as a dark 
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energy particle at late times, accounting for the current phase of cosmic acceleration. We have found that the minimum 
is an attractor with undershoot and overshoot solutions, analogous to normal quintessence models. In studying the 
approach to the attractor, it is important to take into account the kicks due to species becoming non-relativistic 
during the radiation era. These kicks successively push the chameleon field towards the minimum of the potential 
and consequently greatly expand the basin of attraction. It would be interesting to study the role played by the kicks 
in models of quintessence coupled to dark matter. 

The most stringent constraint on our model comes from time variation of masses since BBN. This requires the 
chameleon to have settled to the minimum by the onset of BBN. If this is not realized, then the electron kick will 
result in an unacceptably large variation in masses. This condition is fulfilled for a broad range of initial conditions 
spanning many orders of magnitude in initial chameleon energy density. The allowed range is in fact broader than in 
normal quintessence, largely due to the kicks which make the attractor mechanism comparatively more efficient. 

We have studied the chameleon during inflation and find that it quickly stabilizes at the minimum as the universe 
inflates. We have argued that the chameleon cannot lead to a time variation of the fine-structure constant sufficiently 
large to be consistent with recent observational claims. Finally we applied the chameleon to the cyclic universe. When 
the chameleon mechanism is taken into account the class of potentials relevant to the cyclic universe is enlarged and 
constraints on the parameters relaxed. 
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XI. APPENDIX 
A. The Thin-Shell property 

We have stated in Sec. Ill Al that M < 10^^ eV results from local tests of gravity. We provide more details here 
on how this bound is obtained. Consider a spherical body of homogeneous density pc, radius Rc and total mass 
Mc = 47ri?J?/9c/3, immersed in a homogeneous medium of density poo- We denote by (pc and (poo the field values which 
minimize the effective potential for the respective densities. At short distances, the total force F, gravitational plus 
chameleon-mediated, on a test mass is Q 

F={1 + e)F^ , (89) 

where Fn is Newtonian force and 9 is the fractional force due to the chameleon. For small objects, in a way to be 
made precise below, 9 = 2/3^, the usual answer for a scalar field without potential. For much larger objects, however. 
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one finds that 

9 = 2(3^^2^^^, (90) 

where $c = Mc/SnMp^Rc is the Newtonian potential of the object. Thus, for sufficiently large objects such that 
i4>oc - 4>c)/{'2f3Mpi^c) < 1, one has 9 < 2/3^ and the fifth force is suppressed. If (0oo - (/'c)/(2/3Mpi$c) > 1, on the 
other hand, one simply gets 9 = 2/3^. 

The suppression for large bodies is due to the so-called thin shell effect which can be understood as follows. 
Essentially, only a thin shell under the surface of the object contributes to the cf) pull on a test mass. Indeed, one finds 
that the profile of the chameleon field is such that it is nearly constant up to a radius Rg < Rc, where Rg is given by 

Rc — Rs _ (t>oo — 4>c. /„ X 

Rc ~ 6/3Afpi$e ■ ^ ' 

Whether an object has a thin shell or not, that is whether its chameleon-mediated force is suppressed compared to 
the usual answer, depends on the magnitude of this ratio. The shell is thin if 

Rr — Rs 



< 1 . (92) 



Since 4>c <C 0oo for large density contrast between the object and the ambient matter, this generally reduces to 

< $c . (93) 



Mpi 

For the fiducial potential V{4>) — exp(Af"/(/i") with n ^ C(l), applying this condition for typical test masses used 
in laboratory tests of gravity leads to the constraint 

M < 10"^ eV . (94) 



B. The Instantaneous Kick 



In this Section, we present an analytic approximation to the contribution from each mass threshold to T^' during the 
radiation dominated era, as described in Sec. IIV Al The idea is to approximate each contribution as a delta-function 
source, that is, as an instantaneous kick. Specifically, the equation of motion for is approximated by 

4> + + V^ = -^KHMpidit - to) (95) 

where the right hand side corresponds to a species contributing sharply around time tQ. The constant k is of order 
gi/g^,{mi), in the notation of Sec. IIV XI The order of magnitude of the delta function was estimated using the fact 
that the contribution to is of order {i.e., the r function peaks at a value of order unity) and that the width of 
the r function is of order of a fraction of a Hubble time. 

First, let us consider the regime where is much larger than (pmin- In this case, Eq. ()95|l becomes 

-t- 3H(j) = ~f3KHMpiS{t - to) , (96) 

as the potential is negligible. Assuming the field is initially at rest, the solution is 

t<to: = 0z , (97) 
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and 

t>to: (1 - alf3KHaMpi) [ a^^dt + 0, , (98) 



to 



where ag and Hq are evaluated at Iq. The kick thus results in a jump in from its free evolution of 



A(t> = -all3KHoMpi / a^^di , (99) 

J to 

which, in the radiation-dominated era, converges to 

A(j) = -K/3Afpi . (100) 

Recalling that k ^ gi/g*{iTii), this agrees with Eq. H44|l . 

Next, consider the effect of the kick if the field is close to or at the minimum. In this case, we can linearize the 
effective potential about the minimum and obtain 

4) + 3iJ0 + m^(0 - <^min) = ~(3KHMpid{t - to) ■ (101) 

Letting d(f> = (f> — ^min = a^^/^V': we find 

+ + ^ = -/3Ki?Mpia3/25(i - to) . (102) 

Assuming the field is at the minimum initially, ip — 0, and using 3> H^, the solution reads 

t<to: ^- = 0, (103) 



and 



t>to: ?/> = -- ^-^ sin(m(t - to)) ■ (104) 



Averaging over many oscillations, we find 

Thus the kick makes the chameleon oscillate about the minimum. These oscillations are damped due to the expansion 
of the universe, and the field quickly settles back to the minimum. This is confirmed by solving the equations of 
motion numerically, using the exact form for the r function, as shown in Fig. 

Notice that the amplitude of the oscillations is always small compared to Afpi. In particular, using Eq. (jnH), the 
average variation of masses due to the kick behaves like 



m J / 2m^ \ a 
implying no significant modification to BBN due to the electron kick, for instance. 



<1, (106) 



C. Convergence to the Minimum 

Here we wish to prove Eq. (|55|) directly. If we assume that the scale factor is power-law in time, a{t) ~ t^, then 
Eq. H58fl implies that the mass of the chameleon can be written as 

m{t) = I , (107) 
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le+07 le+08 le+09 le+10 le+11 

1+Z 

FIG. 6: The chameleon starts near the minimum and is subjected to the contribution to from some particle species becoming 
non-relativistic. In this example, the potential is V{4>) ~ /(fp', = 1 and M — \Q MeV, the latter once again chosen to 
be much larger than the required M = 10~^ eV due to numerical limitations. We see that the field oscillates and eventually 
settles back to the minimum. 



where B and q are positive constants. 

Small fluctuations (j)\{t) = (f){t) — (pminit) around the minimum 0min are then governed by 



(Ai it) ~ ( -r^t'-" ] , (109) 



with general solution 

6. (f^ ~ t(i-3p)/2z. ( 

,1-9 

where Zr is a superposition of Bessel functions of order r = (1 — 3p)/2(l — q). Without loss of generality we take 
Zr = Jr- Note that the argument of the Bessel function is large since m{t)t/{l — g) ~ m/H ^ 1. Thus, using the 
asymptotics of the Bessel function, we find 



*.(')~"*-""^/^-(?^). ("0) 



where we have dropped an irrelevant phase factor. 
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Averaging over a few oscillations, while noting that mit) is nearly constant during that time (since rn/m ^ H <^ 
m, as seen from Eq. H57|l). we get 

< >-r3p+'?. (Ill) 

Using Eq. (|107|l . we obtain the desired result: 

■m{t) < <j>l >cx t-^P cx , (112) 
in agreement with Eq. (|55|) . It is straightforward to prove that this also holds for the case q = 1- 
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